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I. INTRODUCTION 

In hadron-hadron collisions at high energies, mas- 
sive lepton pairs are produced via the Drell-Yan pro- 
cess In the standard model, colliding partons from 
the hadrons can interact to form an intermediate W or 
7*/Z vector boson that subsequently decays into a lep- 
ton pair. Initial state quantum chromodynamic (QCD) 
radiation from the colliding partons imparts transverse 
momentum (Pr) to the boson and produces an accom- 
panying final state jet or jets. 

A recent advance in QCD fixed-order perturbative cal- 
culations at 0(a 2 s ) is the evaluation of the Drell-Yan 
cross section that is fully exclusive and differential 
The exclusive cross section includes both the lepton pair 
produced via the W or 7* jZ boson intermediate state, 
and the associated final state partons. It includes finite 
boson widths, boson-lepton spin correlations, and 7 — Z 
interference for the 7* jZ intermediate state. 

The QCD calculation of the Drell-Yan-process cross 
section that is differential in transverse momentum for all 
values of Pr employs a resummation formalism [3| that 
merges fixed-order calculations with an all-orders sum of 
large terms from soft and collinear gluon emissions. The 
dynamics at low Pt is factorized into a calculable pertur- 
bative form factor and a hadron-level, non-perturbative 
one that must be measured. The non-perturbative form 
factor also includes the effect of the intrinsic Pr of par- 
tons in the hadron. Refinement of the phenomenology 
needs precise measurements of the transverse momentum 
differential cross section at low Pr from hadron-hadron 
collisions at various center-of-momentum energies, y/s. 

Previous pp measurements at y/s — 0.63 TeV [4, 5] sup- 
port the resummation formalism, but with limited statis- 
tics. The next pp measurements at y/s = 1.8 TeV 
contributed to the phenomenology at low Pr 0. Re- 
cent pp measurements at y/s — 1.96 TeV [11] are pre- 
cise enough to constrain phenomenological calculations 
of the Drell-Yan lepton pair Pr distribution. Early Large 
Hadron Collider pp results [H, [H| at y/s = 7 TeV show 
agreement with calculations. 

In this article, a new and precise measurement of the 
differential cross section in Pr for Drell-Yan lepton pairs 
from pp collisions at y/s = 1.96 TeV is presented. The 
specific Drell-Yan process is pp — > e + e~ + X, where the 
e + e~ pair is produced through an intermediate 7*/Z bo- 
son, and X is the hadronic final state associated with the 
production of the boson. The measurement of the differ- 
ential cross section is restricted to dielectron pairs within 
the 66-116 GeV/c 2 mass range and is fully corrected to 
include all boson rapidities, electron phase space, and 
detector effects. Within this mass range, the dielectron 
pairs originate mostly from the resonant production and 
decay of Z bosons. 

The cross section, measured using 2.1 fb _1 of collisions 
recorded by the Collider Detector at Fermilab (CDF) 
during 2002-2007, covers < P T < 350 GeV/c. This 
range is subdivided into variable- width Pr bins. For 



P T < 25 GeV/c, the bin width is 0.5 GeV/c. The cross 
section presented for each Pr bin is the average bin cross 
section, Act/ Apr, where Act is the cross section in a Pr 
bin, and APr its width. 

The Act/ APr measurement depends on the correct 
modeling of the physics and detector to unfold the ef- 
fects of the detector acceptance and resolution for the 
pp production of Drell-Yan e + e~ pairs. The modeling 
of the physics and detector is data driven. This mea- 
surement is an extension of the CDF measurements of 
the Droll Yan e + e~ pair rapidity differential cross sec- 
tion [1^1 . and of the decay-electron angular-distribution 
coefficients [lj| that reflect the polarization state of the 
intermediate 7*/^ boson produced in pp —> -f* /Z + X. 
The Act/ APr measurement uses the same 2.1 fb _1 data 
set and analysis methods developed in those measure- 
ments, where both the data and the modeling of the 
physics and detector are well studied and understood. 

Section HT1 provides a brief overview of the QCD calcu- 
lations of Act/ APr used for comparison with this mea- 
surement. Section IlIII provides a summary of CDF and 
the Tevatron collider at Fermi National Accelerator Lab- 
oratory. Section ITVl reports the selection of electrons and 
dielectrons for the Act/ APr measurement. SectionlVlde- 
tails the simulation of the data. Section IVT1 describes the 
cross section and its measurement. Section IVIII is the 
summary. 



II. QCD CALCULATIONS 

For the Drell-Yan process, QCD radiation from the 
colliding partons of the hadrons in the initial state im- 
parts transverse momentum to the lepton pairs. Fixed- 
order perturbative calculations are expected to become 
increasingly reliable with larger transverse momentum. 
However, the Drell-Yan process has two energy scales: 
the lepton-pair invariant mass and transverse momen- 
tum. Difficulties arise in the perturbative calculation 
when these two scales differ significantly. This is a QCD 
multi-scale problem. Simpler perturbative QCD calcula- 
tions usually have one scale, and this scale is often used 
as the scale in the strong coupling, a s , to control accu- 
racy. In addition, all perturbative QCD calculations have 
an arbitrary mass factorization scale that separates the 
hard parton scattering from the soft parton distribution 
functions (PDFs) of the hadrons. With multiple scales, 
scale issues can be harder to control and quantify. 

At the opposite end corresponding to low transverse 
momentum, large contributions from soft and collinear 
gluon emissions begin to dominate and limit the appli- 
cability of standard perturbative calculations. The QCD 
resummation methods are used to overcome this limita- 
tion Q. These resummation methods may be viewed as 
techniques to control large and unreliable contributions 
from multiple QCD scales in the low transverse momen- 
tum kinematic region. 

As neither calculation is expected to be accurate over 
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the entire range of Pr , it is useful to compare them with 
measurements. Of interest is the low Pt region where 
the bulk of events is produced. The understanding and 
proper modeling of QCD at low Pt is important for many 
physics measurements. The Drell-Yan process can be 
used as a benchmark. The measurement presented here 
is compared with a recent QCD resummation calculation, 
RESBOS [nj Hq - [l8l |. and a state-of-the-art QCD fixed- 
order 0(a 2 s ) calculation (NNLO) of A<t/AP t , fewz2 [2, 

The FEWZ2 NNLO calculation is fully exclusive and 
differential for the final-state leptons and partons, and 
includes 7* /Z finite decay width and lepton correlation 
effects. For calculations, the MSTW2008 [H NNLO nu- 
cleon PDFs with their 90% C.L. uncertainties and the 
default FEWZ2 electroweak parameters of the Fermi cou- 
pling constant (G^) scheme and fine-structure constant 
at the Z-boson mass (a~^ = 128) are used. The QCD 
factorization and renormalization scales are both set to 
the Z-boson mass. As no significant phase-space restric- 
tions are applied on the final state, except for the 66- 
116 GeV/c 2 dilepton mass range limit, fewz2 is used 
here as an inclusive calculation. The numerical integra- 
tion accuracy is set to the 1% level. 

The RESBOS calculation utilizes the Collins, Soper, 
and Sterman (CSS) resummation formalism that com- 
bines fixed-order perturbative QCD calculations with an 
all-orders summation of large terms from gluon emis- 
sions [3]. The CSS cross section consists of two terms: a 
W function, which contains the large terms from gluon 
emissions; and a Y function, which is the fixed-order 
cross section minus its asymptotic (large gluon emission) 
terms already in W. The Y function becomes impor- 
tant as the magnitude of the Pr approachs the lepton- 
pair invariant mass. After a Fourier transformation from 
transverse momentum to its conjugate impact-parameter 
space (&), the resummation in the W function is ex- 
pressed as renormalization group equations [2l[. With 
this formalism, the lepton-pair mass and impact param- 
eter scales are connected by the renormalization group 
evolution, through which large perturbative terms are 
reliably controlled. At small b, W is evaluated to ar- 
bitrary order in the renormalized coupling. At large b, 
hadron level, non-perturbative terms that must be mea- 
sured become dominant. The methodologies at small 
and large impact parameters are joined by factorizing W 
into a perturbative and a non-perturbative form factor. 
The perturbative form factor uses the impact parameter, 
6* = b/ y/l + (&/&max) 2 , so that it becomes constant in 
the non-perturbative region. 

The CSS gluon resummation W and Y functions 
should be evaluated to all orders of a s and then com- 
bined to fully describe the physics at all Pt- However, 
practical implementations of the CSS gluon resummation 
formalism evaluate the perturbative Y function and the 
perturbative part of the resummed W function term to 
a finite order in a s . Even with a finite order expansion, 
the CSS gluon resummation formalism provides a good 



description of the physics at low lepton-pair Pp. Above 
a Pp value of about the boson mass, the resummed cross 
section is dominated by the Y function and is close to the 
pure fixed-order calculation. However, in an intermedi- 
ate Pt zone starting from about half the boson mass, 
the cancellation between the W and Y functions eval- 
uated at finite order becomes inadequate because of an 
order mismatch. The W perturbative expansion terms 
are intrinsically all-orders from the underlying resumma- 
tion formalism, but the Y terms are strictly finite-order. 
Within this intermediate Pr zone, W + Y loses accuracy 
and requires compensation in practical implementations 
of the resummation formalism. 

The RESBOS implementations of the W and Y func- 
tions are calculated using CTEQ6.6 PDFs [22|, and are 
provided within RESBOS as cross-section tables on a grid 
of the boson mass, transverse momentum, and rapidity. 
The RESBOS non-perturbative form factor [|(| of the W 
function for the Drell-Yan process is 



exp 



-9i 



1 Q 



gxg 3 ln(100xix 2 ) 



where g x = 0.21 GeV 2 , g 2 = 0.68 GeV 2 , 53 = -0.6, 
Q is the lepton pair mass, Qq = 1.6 GeV/c 2 (with 
frmax — 0.5 GeV -1 for 6* in the perturbative form factor), 
and X1X2 = Q 2 /s. The gi-3 are parameters derived from 
measurements. This form factor describes both low- and 
high-mass data at various y/s from fixed target to collid- 
ers. The specific W and Y function cross-section tables 
used are W321 and Yfc, respectively, and the numerical 
integration uncertainties of RESBOS are under 1% and 
negligible. 

The CSS gluon resummation W function has three sep- 
arate perturbative functions: A, B, and C. In the RES- 
BOS implementation [17| of the W function, W321, those 
functions are evaluated to O(al), C(a 2 ), and 0(a s ), re- 
spectively. Its Y function is 0{a 2 s ). At large Pr, RESBOS 
utilizes both the resummed cross section, W + Y, and 
the C(a 2 ) fixed-order cross section. The resummed cross 
section becomes inaccurate in the intermediate transverse 
momentum region starting from about half of the boson 
mass because of the intrinsic order mismatch described 
previously. Therefore, as the Pr increases, a matching 
procedure between the resummed and fixed-order cross 
section is implemented by RESBOS to provide a reliable 
prediction over all transverse momentum. This match- 
ing is implemented in the Yfc cross-section tabl^E an d is 
a non-trivial, phenomenological part of RESBOS. On the 
other hand, in the transverse momentum region above 
the order of the boson mass, the RESBOS calculation and 
its accuracy are similar to the FEWZ2 NNLO inclusive 



1 To reduce the time needed to compute the cross section 
table to O(a^), the computation is implemented by an 0(a s ) 
calculation with boson mass, rapidity, and transverse momentum 
dependent NNLO-to-NLO K-factors. 
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calculation considered here. The RESBOS calculation also 
includes the full j* /Z interference effects with a finite 
decay-width for the Z boson and with lepton correla- 
tions. The dominant electroweak corrections are included 
in the calculation using the effective Born approximation, 
as done in the LEP electroweak precision measurements. 



III. THE EXPERIMENTAL APPARATUS 

The CDF 11 [23] is a general purpose detector at the 
Fermilab Tevatron Run II pp collider whose center-of- 
momentum energy is 1.96 TeV. The CDF positive z- 
axis is along the proton direction. For particle trajec- 
tories, the polar angle 8 is relative to the proton direc- 
tion and the azimuthal angle <f> is about the beamlinc 
axis. The energy and momentum of a particle are de- 
noted as E and P, respectively. Their components trans- 
verse to the beamline are defined as Et = EsinO and 
Pt = PsinO, repectively. The particle rapidity, y, is 
y = \ ln[ (E + P z c)/(E — P z c)], where P z is the compo- 
nent of momentum along the z-axis. The pseudorapidity 
of a particle trajectory is 77 = — lntan(#/2). Fixed de- 
tector coordinates are specified as (zydet , </>) , where 7yd e t is 
the pseudorapidity from the detector center (z = 0). Por- 
tions of the detector relevant to this analysis are briefly 
described next. 

The central tracker (COT) is a 3.1 m long, open cell 
drift chamber that extends radially from 0.4 m to 1.3 m. 
The 2.1 m long silicon tracker surrounds the Tevatron 
beam pipe and is within the inner radius of the COT. 
Combined, these two trackers provide efficient, high res- 
olution tracking over |?7dct| < 1-3. Both trackers are im- 
mersed in a 1.4 T axial magnetic field produced by a 
superconducting solenoid just beyond the outer radius of 
the COT. 

Outside the solenoid are the central calorimeters, cov- 
ering 1 77det I < 1-1. The forward regions, 1.1 < |?7dot| < 
3.6, are covered by the end-plug calorimeters. All 
calorimeters are scintillator-based sampling calorimeters 
read out with phototubes. Both calorimeters are seg- 
mented along their depth into electromagnetic (EM) and 
hadronic (HAD) sections and transversely into projective 
towers. The EM calorimeter energy resolutions measured 
in test beams with electrons are a/E = 14%/ V-Et for 
the central calorimeter, and a/E = 16%/VP© 1% for 
the plug calorimeter, where the symbol ® is a quadra- 
ture sum, and Et and E are in units of GeV. Both the 
central and plug EM calorimeters have preshower and 
shower-maximum detectors for electromagnetic shower 
identification and shower centroid measurements. The 
combination of the plug shower-maximum detector and 
silicon tracker provides enhanced tracking coverage to 
|?7det|=2.8. 

The Fermilab Tevatron collides bunches of protons and 
anti-protons at a nominal crossing frequency of 2.5 MHz. 
Over 2002-2007 operations, the instantaneous pp colli- 
sion luminosities at the start of collisions increased over 



an order of magnitude to 280 x 10 30 cm _2 s _1 . Colli- 
sion luminosities are continuously measured by the gas 
Cherenkov counters which are just outside the Tevatron 
beam pipe and are in the region 3.7 < |?7dct| < 4.7 [24| . 

The CDF event trigger system has three tiers, LI, L2, 
and L3. The LI trigger is entirely implemented in hard- 
ware, is based on trigger primitives, and is synchronous 
and deadtime-less. Trigger primitives are quantities from 
the front-end readout used for trigger decisions. The L2 
trigger, which processes events selected by the LI trigger, 
is asynchronous and is a combination of hardware and 
software that uses LI primitives along with additional 
front-end data. The L3 trigger processes events selected 
by the L2 trigger and is a speed-optimized version of 
the CDF offline reconstruction. Track- and EM-objects, 
which are available at all trigger levels and are refined at 
each level, form the basis of very efficient trigger paths 
for the electrons used in this measurement. 



IV. DATA SELECTION 

The data set consists of 2.1 fb _1 of pp collisions at 
= 1-96 TeV collected during 2002-2007. Collisions 
producing massive Drell-Yan dielectron pairs have the 
following experimental signatures: 

• A large fraction of the electrons have high Et ■ 

• There are two well-separated electrons of opposite 
charge. 

• The electrons tend to be separated from jets and 
other particles from the interaction. 

These features are used in the selection of events both 
at the trigger and analysis levels. Electrons in both the 
central and plug calorimeters are selected. 

A. Triggers 

The high Et electrons are selected from generic pp 
collisions by two non- attenuated (full-rate) triggers: the 
CENTRAL- 18, and Z-no-track. Each has well-defined 
LI, L2, and L3 trigger paths for both physics and trig- 
ger efficiency measurements. Independent and dedicated 
trigger paths are used for the efficiency measurements. 

The central- 18 trigger is the inclusive electron trig- 
ger for electrons with Ej~ > 18 GeV in the central 
calorimeter region [23| . A track is required at all trigger 
levels. Loose criteria applied at each level select candi- 
dates that are consistent with an electron showering in 
the calorimeter, including EM-shower-like lateral shower 
profile in the EM compartment, EM-shower-like energy 
leakage in the HAD compartment, and matching between 
the track and the shower centroid in the EM shower- 
maximum detector. There is no equivalent inclusive plug 
electron trigger because the LI and L2 tracking and the 
plug calorimeter acceptance do not overlap. 



7 



The Z-NO-TRACK trigger identifies dielectrons using 
solely calorimeter information. No tracking information 
is used. Electron candidates can be in either the central 
or plug calorimeter region. Both candidates are required 
to have Pt > 18 GeV. The only other requirement is 
that shower energy sharing in the EM and HAD com- 
partments be electron-like. While this trigger is specifi- 
cally for dielectron candidates that are both in the plug 
calorimeter region, it accepts the small fraction of dielec- 
tron events that fail the central- 18 trigger. 

B. Electron Selection 

To improve thepurity of the sample, CDF standard 
central and plug [231 ] electron identification requirements 
are applied. Fiducial requirements are always applied to 
ensure that the electrons are in well-instrumented regions 
of CDF where their reconstruction is well understood and 
predictable. Each electron candidate is required to have 
an associated track. Having track matching on both elec- 
tron candidates significantly reduces backgrounds. 

The track vertex position along the beamline (^ v tx) is 
restricted to the inner region of CDF: \Z vtx \ < 60 cm. 
For 2002-2007 Tevatron operations, 4% of the pp lumi- 
nous region along the beamline is outside this fiducial 
region. The pp collision profile along the beamline is 
measured by vertexing multiple tracks in minimum-bias 
events. The multiple track vertexing acceptance is rela- 
tively fiat within |2T v tx| ~ 100 cm - 

As electrons in both the central and plug calorimeter 
regions are used, there are three exclusive Drell-Yan di- 
electron topologies: CC, CP, and PP, where the C (P) 
refers to an electron in the central (plug) calorimeter. In 
the measurement of the ee-pair Pt distribution, the kine- 
matic region of the ee-pair extends over all rapidities, but 
is restricted to the 66-116 GeV/c 2 pair mass range. The 
kinematic and fiducial regions of acceptance for electrons 
in the three dielectron topologies are listed below. 

1. Central-Central (CC) 

• P T > 25 (15) GeV for electron 1 (2) 

• 0.05 < |tjdet| < 1-05 

2. Central-Plug (CP) 

• Et > 20 GeV for both electrons 

• Central region: 0.05 < |?7dct| < 1-05 

• Plug region: 1.2 < |?/ dot | < 2.8 

3. Plug-Plug (PP) 

• Et > 25 GeV for both electrons 

• 1.2 < hdctj < 2.8 

The CC electron Et selection is asymmetric. Electron 
1 has the highest Et- The asymmetric selection is the 
result of an optimization based on the decay electron 



angular distribution measurement [151 ]. It improves the 
acceptance in the electron phase space. The PP electron 
candidates are both required to be in the same end-plug 
calorimeter, and these pairs extend the rapidity coverage 
to \y\ ~ 2.9. At the Tevatron, the kinematic limit for 
\y\ of the ee-pair at the Z-boson mass is 3.1. Drell-Yan 
dielectrons in opposing end plug calorimeters have little 
acceptance, tend to be at low ee-pair rapidities, and are 
overwhelmed by QCD di-jet backgrounds. 

As Drell-Yan high-PT leptons are typically produced 
in isolation, the electron candidates are required to be 
isolated from other calorimetric activity. The isolation 
requirement is that the sum of Et over towers within a 
0.4 isolation cone in (77, <fi) surrounding the electron clus- 
ter be under 4 GeV (E iso < 4 GeV). The towers of the 
electron cluster are not included in the sum. While this is 
a topological selection rather than an electron identifica- 
tion selection, it is included in the electron identification 
efficiencies. 

Electron identification in the central calorimeter region 
is optimized for electrons of Pt > 10 GeV/c. It utilizes 
the COT and silicon trackers, the longitudinal and lat- 
eral (tower) segmentation of the EM and HAD calorime- 
ter compartments, and the shower-maximum strip de- 
tector (CES) within the EM calorimeter. The most dis- 
criminating information is provided by the trackers in 
combination with the CES. An electron candidate must 
have shower clusters within the EM calorimeter towers 
and CES that have EM-like lateral shower profiles. A 
candidate must also have an associated track that ex- 
trapolates to the three-dimensional position of the CES 
shower centroid. The track transverse momentum, Pt, 
must be consistent with the associated electron shower 
Et via an E/P selection when Pt < 50 GeV/c. For 
both the track matching in the CES and E/P selection, 
allowances are included for bremsstrahlung energy loss 
in the tracking volume, which on average is about 20% 
of a radiation length. The fraction of shower energy in 
the HAD calorimeter towers behind the EM tower cluster 
must be consistent with that for electrons (Phad /Pem 
requirement). These selections are more restrictive than 
the ones used in the trigger. 

The central electron selection as described has high pu- 
rity and is called the tight central electron (TCE) selec- 
tion. Its average selection efficiency is 84%. The track- 
finding efficiency on the associated tracks is 99%. To 
improve the selection of central dielectrons, a looser se- 
lection, called the loose central electron (LCE) selection, 
is used on the second electron. The LCE selection does 
not use transverse shower shape constraints, the E/P 
constraint, nor track matching in the CES. For track as- 
sociations, the track need only project into the largest- 
energy calorimeter tower within the cluster of towers as- 
sociated with the EM shower. For electron candidates 
that fail the TCE selection, the LCE selection has an 
average exclusive efficiency of 76%. 

Electron identification in the forward plug calorime- 
ter region also utilizes the COT and silicon trackers, the 
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longitudinal and lateral (tower) segmentation of the EM 
and HAD calorimeter compartments, and the shower- 
maximum strip detector (PES) within the EM calorime- 
ter. However, as the plug calorimeter geometry is com- 
pletely different from the central geometry, the details of 
the identification requirements differ. 

The plate-geometry, end-plug calorimeters have pro- 
jective towers, but these towers are physically much 
smaller than the central calorimetry towers. EM showers 
in the plug calorimeter are clustered into "rectangular" 
3x3 tower clusters in (77, <j>) space, with the highest- 
energy tower in the center. The EM calorimeter energy 
resolution and lateral shower shapes measured in an elec- 
tron test beam use 3x3 shower clustering [25| . The EM 
preshower detector is the first layer of the EM calorime- 
ter and it is instrumented and read out separately. As 
there are ~ 0.7 radiation lengths of material in front of 
it, its energy is always included in the EM-cluster shower 
energy. 

An electron in the plug calorimeter, like those in the 
central region, must also have shower clusters within the 
EM calorimeter towers and PES that have EM-like lateral 
shower profiles. The longitudinal -Ehad/-Eem leakage re- 
quirement is more restrictive because of the deeper depth 
of the EM section and the differing collision conditions in 
the forward region. The plug selection efficiency without 
the tracking requirement averages about 84%. 

Tracks going into the plug calorimeters have limited 
geometrical acceptance in the COT for |r/det| > 1-3. 
The forward tracking coverage of the silicon tracker is 
exploited with a calorimetry-seeded tracking algorithm 
called "Phoenix" . It is similar to central tracking, where 
tracks found in the COT are projected into the silicon 
tracker and hits within a narrow road of the trajectory 
seed silicon track reconstruction. With the Phoenix al- 
gorithm, the track helix in the magnetic field is speci- 
fied by the position of the pp collision vertex, the three- 
dimensional exit position of the electron into the PES, 
and a helix curvature. The curvature is derived from 
the Et of the shower in the EM calorimeter. As the 
Et provides no information on the particle charge, there 
are two potential helices, one for each charge. The al- 
gorithm projects each helix into the silicon tracker and 
seeds the silicon track reconstruction. If both projections 
yield tracks, the higher quality one is selected. The COT 
is not directly used, but tracks found by the trackers are 
used to reconstruct the location of the pp collision vertex. 

The radial extent of the PES, relative to the beamline, 
is 12-129 cm. Depending on the track vertex location 
along the beamline (Z vtx ), a track traverses from to 8 
layers of silicon. A Phoenix track is required to have at 
least three silicon hits. Only plug electrons associated 
to tracks that traversed at least three silicon layers are 
accepted. Eighty percent of the tracks traverse four or 
more silicon layers. Within the plug region, the average 
Phoenix track acceptance is 94% and within this accep- 
tance zone, the track-finding efficiency is 91%. 

The Phoenix algorithm is efficient and results in low 
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FIG. 1. The raw ee-pair Pt distribution for all dielectron 
topologies combined. No corrections or background subtrac- 
tions are applied. The highest Pt is 327 GeV/c. 



background. While the pointing resolution of a Phoenix 
track is good (1 mrad or better), its path length in the 
magnetic field at large |?jdet| is small and the helix curva- 
ture resolution is poor. Consequently, there is neither a 
Pt nor E/P requirement for plug electron identification. 

The central region tracking algorithm utilizes hits 
in the silicon tracker if available. However, the plug 
Phoenix tracking algorithm requires a fully functional 
silicon tracker. This silicon requirement reduces the ef- 
fective integrated luminosity of CP and PP topology di- 
electrons relative to CC dielectrons by 6%. 



C. Dielectron Selection 

Events are required to have a reconstructed dielectron 
pair mass 66 < M ee < 116 GeV/c 2 . For dielectrons of the 
CC topology, the two tracks are required to have opposite 
charge. However, for CP and PP topology dielectrons, 
there is no opposite charge requirement because of the 
significant charge misidentification on Phoenix tracks at 
large |?7 dot |. 

The efficiency for the trigger to select events is typi- 
cally over 99% for dielectrons that pass offline event se- 
lections. The CENTRAL- 1 8 trigger has an inefficiency of 
3% per single central electron due to track association re- 
quirements. The Z-no-track trigger is on average less 
than 0.5% inefficient for all topologies, and complements 
the central- 1 8 trigger. 



D. Measurement Event Sample 

The numbers of events passing all previously described 
selections in the CC, CP, and PP dielectron topologies 
are 51 951, 63 752, and 22 469, respectively. Figure Q] 
shows the raw ee-pair Pt distribution for these events. 
The backgrounds are small, and are from QCD or from 
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WW, WZ, ZZ, tt, W+jets, and Z -> t+t- sources with 
real high-i^x electrons. The QCD background is primar- 
ily from dijets where a track in a jet fakes an electron or 
is an electron from a photon conversion. The high-i?T 
electron sources have at least one real electron. The sec- 
ond electron is either a real second electron or a fake one 
such as in W + jets. 

Overall, the background from QCD and non-Drell-Yan 
high-E'x electrons is 0.5%. It is negligible at low pair 
P T , and for P T > 100 GeV/c, it reaches the 5% level. 
These backgrounds are subtracted from the Pp distribu- 
tion shown in Fig. [T] for the measurement of Aa/APj. 
Backgrounds are significantly reduced, particularly at 
large Pt, by requiring each electron candidate to have 
an associated track. 

The overall QCD background level is 0.3%, and it is 
under 1% at all -Pt- It is estimated with the data used 
for the Pt measurement using an "isolation extrapola- 
tion" procedure. All selection criteria are applied to both 
electron candidates except the isolation energy (E iso ) re- 
quirement on one electron candidate. Its E[ SO distribu- 
tion has a sharp peak at low Ei so from Drell-Yan elec- 
trons (the signal) and a broad, flat distribution extending 
to very large Ei so from QCD sources (the background). 
The E[ SO distribution is fit to a signal plus a background 
component over the full i?i SO range, and the background 
component is extrapolated into the signal region for the 
QCD background estimate. The signal and background 
shapes are derived from the unbiased data set used in the 
measurement, and with selections close to the electron se- 
lections to avoid biases. For the background shape event 
selection, two electron-like candidates are required, but 
one is selected to be "jet-like" by reversing the selection 
requirement on its E- lso and Ehad/Eem parameters. The 
other, whose E lso distribution is the background shape, 
has all electron selection requirements except Ei SO ap- 
plied. 

The high Pt electron backgrounds from WW, WZ, 
ZZ, tt, W + jets, and Z —> t + t~ are derived from the 
simulated samples. The overall background level from 
these sources is 0.2%, but they are the source of the 5% 
backgrounds for Pp > 100 GeV/c. 

Above the Pr of 150 GeV/c, there are 55 events. The 
ee-pair mass distribution has a clear Z-boson mass peak, 
and within the 66-116 GeV/c 2 mass range, there is no 
indication of unexpected backgrounds. The peak location 
and width are consistent with expectations. 



V. DATA SIMULATION 

The acceptance for Drell-Yan dilepton pairs is ob- 
tained using the Monte Carlo physics event generator, 
PYTHIA 6.214 and the CDF event and detector sim- 
ulations, pythia generates the hard, leading order (LO) 
QCD interaction, q + q — > 7* /Z, simulates initial state 
QCD radiation via its parton shower algorithms, and gen- 
erates the decay, j*/Z ->■ The CTEQ5L [27| nu- 



cleon parton distribution functions (PDFs) are used in 
the QCD calculations. The underlying event and "/* /Z 
boson Pp parameters are pythia tune AW (i.e., PYTUNE 
101, which is a tuning to previous CDF data) |26l. |28|. 

Generated events are processed by the CDF event and 
detector simulation. The event simulation includes PHO- 
TOS 2.0 [2{| which adds final-state quantum electrody- 
namics (QED) radiation to decay vertices with charged 
particles, e.g. 7* /Z —> e + e~ . The time-dependent beam 
and detector conditions for data runs recorded and used 
for physics analyses are simulated. The beam conditions 
simulated are the p and p beamline parameters, the pp 
luminous region profile, and the instantaneous and in- 
tegrated luminosities per run. The detector conditions 
simulated are detector component calibrations, which in- 
clude channel gains and malfunctions. Thus, the simu- 
lated events parallel the recorded data, and are recon- 
structed, selected, and analyzed as the data. 

The Aer/APr measurement is data driven and depends 
on the correct modeling of both the physics and the de- 
tector. The procedure involves the measurement and 
tuning of the underlying kinematics and detector param- 
eters that make the simulated, reconstructed event dis- 
tributions match the actual data as precisely as possible. 
This is a bootstrap process that iterates if necessary for 
the required precision. The default simulation does not 
reproduce the data at the precision required. The follow- 
ing subsections describe the model tunings. 



A. Physics Simulation 

The Drell-Yan dilepton production is described by 

d 4 a _ d 3 a dN 

dM 2 dy dP T dQ ~ dM 2 dy dP T dSl ' 

where d 3 a/dM 2 dy dPj- is the unpolarized j*/Z boson 
production cross section at the resonance mass M with 
subsequent decay to e + e~, and dN/d£l the electron an- 
gular distribution of the 7* /Z — > e + e~ decay in a boson 
rest frame. For this measurement, the single differential 
distributions da / dy and da / dPr , and the electron angu- 
lar distribution, are tuned to the data. The y distribution 
tuning for ^* / Z production is from the da/dy measure- 
ment [l4j . The tuning of the electron angular distribution 
is briefly reviewed next. The tuning of da/dPr is specific 
to this analysis, and is presented last. 

The pythia parton showering starts with the qq — > 
7*/Z annihilation vertex at the end of the shower chain 
then evolves the shower backwards in time to an ini- 
tiating qq or qg state. The Compton production pro- 
cess cannot be fully simulated. While its gluon split- 
ting subprocess is simulated, the gluon fusion subprocess, 
qg —>g*—> 9 + 7* /Z, cannot be simulated from the anni- 
hilation vertex at the end of the shower chain. The gluon 
fusion production rate is compensated in the shower, but 
there is no compensation to the boson polarization states 
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affected by this subprocess. The boson polarization af- 
fects the decay electron angular distribution. 

The decay electron angular distribution is analyzed in 
the Collins-Soper (CS) rest frame [30| of the e + e~ pair. 
The CS frame is reached from the laboratory frame via a 
Lorentz boost along the lab z-axis into a frame where the 
z-component of the pair momentum is zero, followed by 
a boost along the Pr of the pair. At Pr = 0, the CS and 
laboratory coordinate frames are the same. Within the 
CS frame, the z-axis for the polar angle is the angular 
bisector between the proton direction and the negative 
of the anti-proton direction. The x-axis is the direction 
of the Pt- The polar and azimuthal angles of the e~ in 
the rest frame are denoted as $ and <p, respectively. 

The general structure of the Drell-Yan decay lepton 
angular distribution in a boson rest frame consists of nine 
helicity cross sections governed by the polarization state 
of the vector boson 13111 . 



16tt dN 



= (1 + cos 2 I?) + 

A i (1 - 3cos 2 tf) + 

Ai sin I'd cos ip + 

A2 — sin 2 d cos 2(p + 

A3 sin d cos tp + 

A4 cos d + 

A§ sin 2 ■& sin 2(p + 

Aq sin I'd sin (p + 

A-j sin i? sin ip . 



The Aq_7 coefficients are cross section ratios, and are 
functions of the boson kinematics. They are zero at 
Pr = 0, except for the electroweak part of A4 respon- 
sible for the forward-backward e~ asymmetry in cos$. 
The ^5_7 coefficients appear at C(a 2 ) and are small in 
the CS frame. The decay-electron angular-distribution 
analysis [l5| in the CS frame measures the large and ac- 
cessible decay electron angular coefficients, Aq, A2, A3, 
and A4,, as functions of Pt- These measurements are 
incorporated into the modeling of 7* JZ — > e + e~ decays. 

The generator-level Pt distribution is adjusted, bin by 
bin, so that the shape of the reconstruction-level, sim- 
ulated Pr distribution is the same as in the data. The 
method uses the data-to-simulation ratio of the number 
of reconstructed events in Pp bins as an iterative adjust- 
ment estimator for the generator level Pr bins. Succes- 
sive iterations unfold the smearing of events across Pr 
bins. Figure [5] is the generator-level Pr correction func- 
tion that makes the data-to-simulation ratio uniform. 
Statistical fluctuations in the ratio are smoothed out. 
The AiV/APr correction is the measurement of the shape 
of da 1 'dP? that is used in the physics model. 
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FIG. 2. The Pt correction function applied to the generator 
level AN/APt distribution that makes flat the ratio of the 
observed data to the simulated data. The points are at the 
center of the Pt bins. For the low-statistics Pt > 120 GeV/c 
region, an average correction is used. 



B. Detector Simulation 

The simulation is used to calculate the combined de- 
tector acceptance (A) and selection efficiency (e) as a 
function of kinematic variables for Drell-Yan dielectrons. 
The combined acceptance and efficiency convolution is 
denoted as A<E>e. Single-electron selection efficiencies are 
measured and incorporated into the simulation as event- 
weight scale factors. The scale factors are ratios of the 
measured efficiencies of the data to the simulated data. 

The electron-trigger efficiencies have an Pt (calorime- 
try) and ?7dct (tracking) dependence that are measured 
and incorporated into the simulation. The electron- 
identification efficiencies are measured as a function of 
?7dot for both central and plug region electrons. Plug re- 
gion efficiencies are measured separately for CP and PP 
topology dielectrons due to their different environments. 
Plug-electron efficiencies have a clear time dependence 
due to the increasing instantantaneous luminosities de- 
livered by the Tevatron. This dependence is incorporated 
into the simulation. Luminosity effects are measured us- 
ing the number of pp vertices reconstructed by the track- 
ers per event. 

A precise model of the calorimeter response in the sim- 
ulation is important for the calculation of A®t. Electron 
kinematics are derived from a three-momentum that uses 
the electron energy measured in the calorimeters for the 
momentum magnitude and the associated track for the 
direction. The simulated electron energy scale calibra- 
tion and resolution versus 77det are tuned using the elec- 
tron Pt distribution. The default scale and resolution 
per rjdet bin are adjusted so that the electron Pt distri- 
bution reconstructed in simulation matches that of the 
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FIG. 3. The overall CC topology central electron Et distri- 
bution. The crosses are the data and the histogram is the 
simulated data. 
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FIG. 4. The overall CP-topology plug electron Et distri- 
bution. The crosses are the data and the histogram is the 
simulated data. 



data. Only the constant term in the energy resolution is 
adjusted. Since the default simulation parametrization of 
the energy resolution can already have a constant term, 
the resolution adjustment is done with an additional con- 
stant term C2, 



3500 



a 

E 



' Ci © c 2 



where a is the energy resolution, E is the energy, Cq and 
ci are the default parameters of Section and the © 
denotes combination in quadrature. The tuned values of 
c 2 on average are 0.9% and 2.3% for the central and plug 
calorimeters, respectively. The steeply rising and falling 
parts of the electron Et distribution dominate the con- 
straints. The three dielectron topologies, CC, CP, and 
PP, provide multiple and independent central and plug 
electron Et samples. The r^det-dependent Et distribu- 
tions of each topology are calibrated independently. Af- 
ter the 77det-dependent parameters are determined, the 
separate CC, CP, and PP dielectron mass distributions 
are used to set an overall global scale and resolution ad- 
justment for central and plug electrons. 

The simulation is compared to data using his- 
togrammed electron Et and ee-pair mass distributions. 
Since the backgrounds are small, they are ignored. The 
comparison statistic is the x 2 between the simulation and 
data. The event count of the simulated data is normal- 
ized to that of the data, and only statistical uncertainties 
are used in the calculation. 

The rydct-dependent calorimeter response tunings pro- 
vide a good match between the simulated-data and data. 
Figure [3] shows the Et distribution of CC-topology cen- 
tral electrons. The corresponding plot for PP-topology 
plug electrons is similar in shape except that the width 
of the Et "peak" is slightly narrower. Figure 2] shows 
the Et distribution of CP-topology plug electrons. The 
corresponding plot for CP-topology central electrons is 
very similar. A % 2 test is used to evaluate the compati- 
bility between the simulation and data. For CC-central, 




100 105 110 115 



CC ee-Pair Mass (GeV/c 2 ) 

FIG. 5. The overall CC topology ee-pair mass distribution. 
The crosses are the data and the histogram is the simulated 
data. 



CP-central, CP-plug, and PP-plug electrons, the x 2 val- 
ues are 117, 100, 87, and 135, respectively, for 100 bins 
(90 bins for PP). 

Figure [5] shows the CC-topology ee-pair mass distri- 
bution. The ee-pair mass distributions for the CP and 
PP topologies are similar. The simulated-data to data 
X 2 for the CC-, CP-, and PP-topology ee-pair mass dis- 
tributions are 107, 123, and 114, respectively, for 100 
bins. The sharp and narrow Z-peaks provide significant 
constraints on the the global energy scale and resolution 
parameters. 



VI. THE CROSS SECTION 

The differential cross section in Pt is the average cross 
section in a Pt bin, or Acr/APr, where Act is the inte- 
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grated cross section in a bin. The Aer is defined as: 



G C A<®e 

where N is the background subtracted event count, C 
is the effective integrated luminosity, and A eg e is the 
combined acceptance and efficiency. The effective lumi- 
nosity, C, is 2057 pb^ 1 , and it includes the acceptance of 
the |Z vtx | < 60 cm fiducial restriction. The details of the 
measurement and its uncertainties are presented next. 



A. Acceptance and Efficiency Unfolding 

The combined acceptance and efficiency, A eg e, is 
calculated using the simulation to convolve individual 
electron rj^et acceptances and efficiencies into an ee-pair 
Pr quantity. The value of A eg e ranges from 0.22 at 
P T ~ 0.2 GeV/c to 0.30 at P T ~ 200 GeV/c. As P T 
increases, the ee-pair rapidity becomes more central, the 
electron Pr becomes larger, and the acceptance slowly 
increases. 

The smearing of the observed Pt away from the gen- 
erator (7* /Z) level value is significant relative to the bin 
size at low Pt: It has an rms width of about 2.2 GeV/c 
and is non-Gaussian. Detector resolution and QED radi- 
ation from the J*/Z — > e + e~ vertex induce distortions to 
the reconstructed ee-pair mass and Pt distributions. In 
addition, they induce a broad enhancement in the A (g e 
function. It rises from 0.22 at Pt ~ 0.2 GeV/c to a 
broad maximum of 0.28 around a Pt of 8 GeV/c, then 
decreases to 0.24 at Pr ~ 30 GeV/c before increasing 
again at larger Pr due to the increased acceptance. 

When A <g e is used to calculate cross sections, it un- 
folds the effects of smearing. The (A (g e) _1 correction 
is applied bin-by-bin and consists logically of two steps. 
The first step is a scaling correction on the number of 
reconstructed and selected events. This scales (unfolds) 
the number of events reconstructed in a Pr bin into the 
number of reconstructed events produced in the bin. The 
simulation provides an average scaling factor. The sec- 
ond step is a standard detector acceptance correction on 
this scaled (unfolded) event count. 

For the cross section uncertainty evaluation, more in- 
formation on event production and migration among the 
Pr bins is required. The number of events produced in 
each bin have statistical fluctuations. With smearing, 
there is event migration among the bins, and this mi- 
gration is also subject to statistical fluctuations. At low 
Pr, event migration between bins is large. If these mi- 
grations are unaccounted, the cross section uncertainty 
will be significantly underestimated. As the event migra- 
tion between bins is not measured, these migrations are 
estimated with the simulation. 



B. Unfolding Uncertainty Model 

Comparisons of fully-corrected cross section measure- 
ments with theoretical cross sections are not straightfor- 
ward. Where detector smearing is significant, there are 
significant uncertainty correlations among the Pr bins 
due to the event migrations among the bins. The sim- 
ulation behind the scaling correction accounts for these 
migrations. The scaling correction uncertainty has both 
statistical and systematic biases. The systematic bias is 
from the residual simulation model bias on the bin scal- 
ing factor. This bias has been mitigated by the model 
tuning described in Section IV Bl The sources of statisti- 
cal uncertainty from event migration for the scaling cor- 
rection are discussed, and a model of per-measurcmcnt 
(per-single-experiment) fluctuations for the uncertainty 
that uses the simulation is specified. 

Within the context of the simulation, information 
about the event migration of reconstructed events among 
Pr bins is in its transfer matrix, fiik, where fiik is the ex- 
pectation value of the number of events produced in bin k 
that migrate into bin I. The expectation value of a quan- 
tity is denoted with an overbar, e.g. n. All expectation 
values are normalized to the integrated luminosity of the 
data. The number of events that do not migrate out of a 
bin is denoted by n g . The number of events that migrate 
out and in are denoted by n Q and fii, respectively. In 
terms of the transfer matrix, n/fe, the n g , n , and fii for 
Pr bin m are, respectively, n mm , the sum of ni m over the 
migration index I excluding bin to, and the sum of n m k 
over the production index k excluding bin m. The per- 
measurement statistical fluctuation of a quantity from its 
expectation value is denoted by S followed by the quan- 
tity, e.g., Sn = n — n. An ensemble variance is denoted 
by 6 2 , e.g., for Poisson statistics, 5 2 n = n, and if c is a 
constant, 5 2 cn = c 2 n. 

The scaling correction factor is p = N g /N r , where 
N g = n g + fi is the expectation on the number of events 
produced in a bin, and N r = n g + fii is the number 
of events reconstructed in a bin. Any residual model 
systematic bias is in p. For a given measurement, the 
number of events produced and reconstructed in a bin 
are N g = n g + n a and N r = n g + m, respectively. The 
scaling correction estimate for N g is p N r . The difference 
between the scaling correction estimator p N r and its tar- 
get N g gives a bias between them, B = p N r —N g . If there 
are no target fluctuations (N g = N g ), B is the statistical 
bias of the estimator. With target fluctuations, there are 
two statistical biases, pN r — pN r (— 6p N r ) and N g — N g 
(= SN g ), and B is their difference. 

With no smearing, the estimator and target, along with 
their fluctuations, are identical, so B = and the statis- 
tical uncertainty of the scaling correction is just that of 
the estimator. With smearing, the estimator and target 
fluctuations are not fully correlated, so B ^ and the 
scaling correction statistical uncertainty is from a com- 
bination of estimator and target statistical fluctuations. 
The estimator (pN r ) and target (N g ) have three sta- 



13 













R q 


































- 












































































- 
















































P 






















































































, 1 







e + e' Pair P T (GeV/c) 

FIG. 6. Rg(Pr) and p(Pi). The solid histogram is R g , the 
bin variance of the uncertainty model relative to the variance 
of the produced events. The abrupt drops are where the bin 
size changes. The lower, dashed histogram is p. 

tistically independent elements: n g , n , and rii. As n g 
is part of both the estimator and target, the common 
overlap must be removed to avoid double counting. The 
total per-measurement fluctuation for the scaling correc- 
tion, denoted as SN' g , is defined as the sum of fluctuations 
(8n — n — n) from the estimator and the target minus 
their common term, Sn g : 

5N' g =5pN r + 5N g -Sn g 

— 8 p(n g + m) + S(n g + n a ) — Sn g 
= S p(n g + m) + Sn . 

For Poisson statistics, the Pr bin ensemble variance is: 

8 2 N' g = p 2 (n g +n l )+n 
= pN g + n Q , 

where N g = pN r = p(n g + Hi) is used. The covariance 
from the Sn a and 5 prii terms between bins k and / is 
Pkfiki + pmik- 

The ratio, R g , of S 2 N' g to 5 2 N g (= N g ) is the variance 
of the model relative to the variance of only the produced 
events. Figure [6] shows both the ratio and the scaling 
correction factor as functions of Pt- In the low Pp bins, 
n a and n, are separately much larger than n g . Their 
effects are significant as R g — p + n /N g . 

For the uncertainty evaluations, the cross section is 
rewritten as Act = pN r /(L A'), where A' = p A ® e. 
The uncertainty on C is systematic and is considered 
separately. Thus, the fractional uncertainty on Act is 
a combination of the fractional uncertainty of pN r and 
A'. The fractional uncertainty of pN r is defined as the 
uncertainty of pN r from the model (SN' ) divided by pN r 
(= N g ). The correlation of these fractional uncertainties 
between Pp bins I and k is given by the fractional covari- 
ance matrix: Vik/(N g i N g k), where Vik is the covariance 
matrix of the model, and N g i and N g k are the N g of bin I 



and fc, respectively. The small acceptance fractional un- 
certainties are added in quadrature to the diagonal part 
of the fractional covariance matrix. The measured cross 
sections are used to convert the unitless fractional ma- 
trix into units of cross section squared, and this matrix 
is used to propagate uncertainties for the total cross sec- 
tion measurement and for the comparison of a prediction 
with the measured cross section. 



C. Systematic Uncertainties 

The largest source of uncertainty is the effective in- 
tegrated luminosity, C. It has an overall uncertainty of 
5.8% that consists of a 4% uncertainty of the acceptance 
of the gas Cherenkov luminosity detector [24j to pp in- 
elastic collisions and a 4.2% measurement uncertainty. It 
is common to all Pr bins and not explicitly included. The 
acceptance uncertainty is primarily from the uncertainty 
in the beamline and detector geometry (material), and 
from the uncertainty in the model of the inelastic cross 
section. The inelastic cross section model contributes 2% 
to the acceptance uncertainty. The measurement uncer- 
tainty contains the uncertainty of the absolute pp inelas- 
tic cross section. 

The uncertainty on i®e has a component from the 
input electron efficiency measurements which depend on 
?7dct and instantaneous luminosity. The simulation is 
used to propagate these electron measurement uncertain- 
ties into an uncertainty for the ee-pair Pp and to include 
correlations of the same measurements. The calculated 
uncertainty is uniform and amounts to about 1% over 
< Pr < 20 GeV/c. It slowly decreases at higher P T . A 
large fraction of the uncertainty is due to plug electron 
measurement uncertainties. The fractional uncertainty 
decreases with Pp because the fraction of plug events 
decreases. Because the same measurements are used on 
all Pr bins, the uncertainty is treated as fully correlated 
across bins. 

The calorimeter response modeling uncertainty analy- 
sis is limited by the statistical precision of the simulated 
data. At the peak of the Pr distribution, the statistical 
uncertainty is 0.3%. The variations on the central and 
plug calorimeter global energy scale and resolutions tun- 
ings allowed by the data propagate into changes of A ® e 
that are no larger than its statistical uncertainty. These 
changes are not independent. 



D. Results 

The Droll- Yan Act/APt for e + e~ pairs in the Z-boson 
mass region of 66 — 116 GeV/c 2 is shown in Fig. [7] and 
tabulated in Table fl] The total cross section from the nu- 
merical integration of the cross section in each Pr bin is 
256.1 ± 1.3 ± 2.6 pb, where the first uncertainty is statis- 
tical and the second is the systematic uncertainty due to 
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e + e Pair P T (GeV/c) 

FIG. 7. The Atr/APr cross section versus Pr. Cross section values are plotted at the bin center. The horizontal bars 
represent the bin extent and the vertical bars are the cross section uncertainties. The solid (black) crosses are the data and all 
uncertainties except the integrated luminosity uncertainty are combined and plotted. The solid (red) histogram is the RESBOS 
calculation. The dot-dash (blue) bars of the Pt > 25 GeV/c region are the fewz2 calculation. For the calculations, only 
numerical uncertainties are included but they are too small to be visible. The inset is the Pr < 25 GeV/c region with a linear 
ordinate scale. 



electron efficiency measurements. The 5.8% integrated 
luminosity uncertainty of 14.9 pb is not included. 

Figure [7] shows that the RESBOS prediction has a gen- 
eral agreement with the data over the full range of Pt. 
The RESBOS total cross section from the numerical in- 
tegration of its cross section in each Pt bin is 254 pb. 
Figure [5] shows the ratio of the measured cross section to 
the RESBOS prediction in the lower Pt region. 

The detector smearing correlates neighboring Pr bin 
uncertainties that are estimated with the model specified 
in Section IVI Al For the low Pt bins, the correlations 
spread across many bins but for Pr > 40 GeV/c, the 
correlations are predominantly between nearest neigh- 
bors. The cross section covariance matrix eigenvalues 
and eigenvectors are used for the y 2 comparison between 
the data and RESBOS. The eigenvalues are the measure- 
ment uncertainties (variances) of the associated eigen- 
vector. Measurement uncertainties between eigenvectors 
are uncorrelated. As an eigenvector corresponds to many 
Pp bins because of smearing, its most probable Pr bin 
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FIG. 8. The ratio of the measured cross section to the res- 
BOS prediction in the Pr < 130 GeV/c region. The RESBOS 
total cross section is normalized to the data. The inset is an 
expansion of the low Pr region. 
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TABLE I. The Act/APt cross section versus Pp. The first uncertainty is the statistical uncertainty. The second uncertainty is 
the efficiency measurement systematic uncertainty, which is 100% correlated across all bins. The 5.8% luminosity uncertainty 
applies to all bins but is not included. 



Pt bin 
GeV/c 


Acr/APr 
pb / GeV/c 






Pt bin 
GeV/c 


Acr/APr 
pb / GeV/c 




0.0-0.5 


(3.613 ±0.168 ±0.035) 


X 


10 u 


20.5-21.0 


(2.923 ±0.143 ±0.030 


x 10 u 


0.5-1.0 


(1.008 ± 0.027 ± 0.010) 


X 


10 1 


21.0-21.5 


(2.877 ±0.144 ±0.030 


x 10° 


1.0-1.5 


(1.551 ± 0.033 ±0.015) 


X 


10 1 


21.5-22.0 


(2.603 ±0.134 ±0.027 


x 10° 


1.5-2.0 


(1.947 ±0.037 ±0.019) 


X 


10 1 


22.0-22.5 


(2.624 ±0.139 ±0.027 


x 10° 


2.0-2.5 


(2.158 ± 0.039 ±0.021) 


X 


If) 1 


22.5-23.0 


(2.590 ±0.140 ±0.026 


x 10° 


2.5-3.0 


(2.295 ± 0.040 ± 0.023) 


X 


If) 1 


23.0-23.5 


(2.516 ±0.139 ±0.026 


x 10° 


3.0-3.5 


(2.258 ± 0.039 ±0.022) 


X 


10 1 


23.5-24.0 


(2.200 ±0.124 ±0.022 


x 10° 


3.5-4.0 


(2.235 ± 0.039 ± 0.022) 


X 


10 1 


24.0-24.5 


(1.948 ±0.113 ±0.020 


x 10° 


4.0-4.5 


(2.061 ± 0.037 ±0.021) 


X 


If) 1 


24.5-25.0 


(2.179 ±0.129 ±0.022 


x 10° 


4.5-5.0 


(1.987 ± 0.036 ± 0.020) 


X 


If) 1 


25.0-26.0 


(2.032 ±0.085 ±0.021 


x 10° 


5.0-5.5 


(1.876 ± 0.035 ±0.019) 


X 


10 1 


26.0-27.0 


(1.736 ±0.076 ±0.018 


x 10° 


5.5-6.0 


(1.729 ±0.034 ±0.017) 


X 


10 1 


27.0-28.0 


(1.633 ±0.075 ±0.016 


x 10° 


6.0-6.5 


(1.563 ±0.032 ±0.016) 


X 


If) 1 


28.0-29.0 


(1.616 ±0.077 ±0.016 


x 10° 


6.5-7.0 


(1.468 ±0.031 ±0.015) 


X 


If) 1 


29.0-30.0 


(1.381 ±0.069 ±0.014 


x 10° 


7.0-7.5 


(1.374 ± 0.030 ±0.014) 


X 


10 1 


30.0-32.0 


(1.284 ±0.045 ±0.013 


x 10° 


7.5-8.0 


(1.307 ±0.030 ±0.013) 


X 


10 1 


32.0-34.0 


(1.005 ±0.039 ±0.010 


x 10° 


8.0-8.5 


(1.183 ± 0.028 ±0.012) 


X 


If) 1 


34.0-36.0 


(8.769 ± 0.361 ± 0.088) 


x 10 _1 


8.5-9.0 


(1.112 ±0.027 ±0.011) 


X 


If) 1 


36.0-38.0 


(7.959 ± 0.352 ± 0.079) 


x 10 _1 


9.0-9.5 


(1.033 ± 0.026 ±0.011) 


X 


10 1 


38.0-40.0 


(7.068 ± 0.336 ± 0.070) 


x 10 _1 


9.5-10.0 


(1.024 ±0.027 ±0.011) 


X 


10 1 


40.0-44.0 


(5.605 ±0.193 ±0.055) 


x 10 _1 


10.0-10.5 


(9.043 ± 0.244 ± 0.094) 


X 


10" 


44.0-48.0 


(4.600 ±0.179 ±0.044) 


x 10 _1 


10.5-11.0 


(8.295 ± 0.231 ± 0.084) 


X 


10" 


48.0-52.0 


(3.552 ±0.156 ±0.033) 


x 10 _1 


11.0-11.5 


(8.319 ± 0.239 ±0.085) 


X 


10° 


52.0-56.0 


(2.760 ±0.136 ±0.025) 


x 10 _1 


11.5-12.0 


(7.780 ± 0.229 ± 0.079) 


X 


10° 


56.0-60.0 


(2.311 ±0.128 ± 0.020) 


x 10 _1 


12.0-12.5 


(7.465 ± 0.227 ± 0.076) 


X 


10° 


60.0-65.0 


(1.618 ±0.089 ±0.014) 


x 10 _1 


12.5-13.0 


(6.839 ±0.215 ±0.069) 


X 


10° 


65.0-70.0 


(1.343 ±0.084 ±0.011) 


x 10 _1 


13.0-13.5 


(6.411 ± 0.208 ±0.065) 


X 


10° 


70.0-75.0 


(1.094 ±0.078 ± 0.009) 


x 10 _1 


13.5-14.0 


(6.220 ± 0.208 ± 0.064) 


X 


10° 


75.0-80.0 


(8.415 ± 0.678 ± 0.068) 


x 10~ 2 


14.0-14.5 


(5.890 ± 0.204 ± 0.060) 


X 


10" 


80.0-85.0 


(6.347 ± 0.565 ± 0.049) 


x 10~ 2 


14.5-15.0 


(5.363 ±0.190 ±0.055) 


X 


10" 


85.0-90.0 


(4.982 ± 0.504 ± 0.038) 


x 10 -2 


15.0-15.5 


(5.186 ±0.190 ±0.053) 


X 


10° 


90.0-95.0 


(3.786 ± 0.422 ± 0.028) 


x 10~ 2 


15.5-16.0 


(4.792 ±0.181 ±0.049) 


X 


10° 


95.0-100.0 


(2.988 ± 0.389 ± 0.023) 


x 10 -2 


16.0-16.5 


(4.431 ±0.172 ±0.045) 


X 


10" 


100.0-110.0 


(2.298 ±0.227 ±0.016) 


x 10 -2 


16.5-17.0 


(4.149 ±0.165 ±0.042) 


X 


10" 


110.0-120.0 


(1.449 ±0.178 ±0.010) 


x 10 -2 


17.0-17.5 


(4.346 ±0.179 ±0.044) 


X 


10° 


120.0-130.0 


(9.369 ± 1.389 ± 0.064) 


x 10 -3 


17.5-18.0 


(3.931 ±0.166 ±0.040) 


X 


10° 


130.0-140.0 


(8.395 ± 1.496 ± 0.055) 


x 10 -3 


18.0-18.5 


(3.757 ±0.163 ±0.038) 


X 


10" 


140.0-150.0 


(5.304 ±1.174 ±0.034) 


x 10~ 3 


18.5-19.0 


(3.753 ±0.167 ±0.038) 


X 


10° 


150.0-175.0 


(1.861 ±0.331 ±0.012) 


x 10 -3 


19.0-19.5 


(3.586 ±0.163 ±0.036) 


X 


10° 


175.0-200.0 


(5.283 ± 1.478 ±0.031) 


x 10~ 4 


19.5-20.0 


(3.303 ±0.154 ±0.034) 


X 


10° 


200.0-250.0 


(2.838 ± 1.019 ± 0.019) 


x 10~ 4 


20.0-20.5 


(2.952 ±0.142 ±0.030) 


X 


10° 


250.0-350.0 


(1.489 ±1.162 ±0.009) 


x 10- 4 



is used for its association to a Pt bin. The mapping of 
eigenvectors to Pt bins is described next. 

The Pt bins are numbered consecutively, 0-81 (low- 
est to highest Pt), and the bin number is denoted by n. 
The bin-number expectation values of the eigenvectors 
are used for their assignment to Pr bins. The eigen- 
vector with the lowest expectation value is assigned to 
Pr bin 0, the next lowest to Pr bin 1, and so on. For 
Pt < 25 GeV/c, the rms width of the expectation value 
is about 4 bins, and above it, about 1 bin or less. In the 
13-18 GeV/c region, the rms width is the largest, 5-6 



bins. 

The x 2 is calculated for the eigenvector associated with 
the Pr bin n. For reference, the uncorrelated \ 2 is also 
calculated for the bin. The cumulative x 2 from bin to n 
inclusive is denoted as x 2 { n )- The number of degrees of 
freedom of x 2 ( n ) is n - A useful measure is x 2 ( n )~ n: it i s 
typically constant when the prediction is compatible with 
the data and increases over regions with discrepancies. 

Figure [5] shows the x 2 (n) — n of the the RESBOS predic- 
ton. For the correlated x 2 , changes in x 2 ( n ) ~~ n can only 
be associated with a Pr region because of smearing. In 
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from the measurement is significant, the difference pro- 
vides a measure of the importance of higher order contri- 
butions above O{ot 2 s ). The PDF uncertainties provided 
by FEWZ2 are at the 3% to 4% level. The uncertainties 
from variations of the QCD factorization and renormal- 
ization scales (from the Z-boson mass) in the Pr regions 
of 25-30, 100-110, and 200-250 GeV/c, are at the 7%, 
5%, and 6% level, respectively. However, the accuracy 
of these scale uncertainties is unclear because of the two 
different scales (lepton-pair mass and transverse momen- 
tum) inherent in this QCD calculation. 



VII. SUMMARY 



FIG. 9. X 2 (n) 

— n versus Pt bin number of the RESBOS predic- 
tion. The solid (red) histogram includes bin correlations and 
the dashed (blue) histogram does not. Bins 0-49 cover the 
0-25 GeV/c region, bin 60 is 40 GeV/c, bin 70 is 85 GeV/c, 
and bin 80 is 200 GeV/c. 
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The transverse momentum cross section of e + e~ pairs 
in the Z-boson mass region of 66-116 GeV/c 2 produced 
in pp collisions is measured using 2.1 fb _1 of Run II data 
collected by the Collider Detector at Fermilab. The mea- 
surement is data driven and corrected for the detector ac- 
ceptance and smearing. The physics and detector models 
of the simulation used for the correction are tuned so that 
the simulation matches the data. The precision of the 
data and the measurement method require both the data 
and simulation to be well calibrated and understood. The 
measurement uncertainties are from a simulation-based 
model that quantifies the effects of event migration be- 
tween measurement bins due to detector smearing. 

Comparisons of this measurement with current quan- 
tum chromodynamic 0(a 2 s ) perturbative and all-orders 
gluon resummation calculations show reasonable agree- 
ment. The data is of sufficient precision for further re- 
finements in the phenomenology of the Drell-Yan lepton 
pair transverse momentum distribution. 



FIG. 10. The ratio of the measured cross section to the fewz2 
prediction in the 25 < Pt < 250 GeV/c region. The fewz2 
cross section is not normalized to the data. 



the Pt < 25 GeV/c region (bins 0-49), there are small 
differences but the data may allow further tuning of the 
RESBOS non-perturbative form factor that is important 
in this region. 

In the 44 < P T < 90 GeV/c region of Fig. E (bins 
61-70 of Fig. |9]), the RESBOS prediction is systematically 
lower than the data. This region is where the resummed 
calculation must be matched to the fixed-order perturba- 
tive calculation. This region is where the data can also 
contribute to the RESBOS resummation phenomenology 
of the Drell-Yan lepton pair Pt distribution at the Teva- 
tron. 

Figure [TU] shows the ratio of the measured cross section 
to the FEWZ2 prediction. There is reasonable agreement 
with the data in the high Pp region where the RESBOS 
and fewz2 calculations are in agreement with each other. 
In Pt bins where the deviation of the fewz2 prediction 
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